A new mechanism for generating broadband 
puisar-iilce poiarization 



Houshang Ardavan/ Arzhang Ardavan,^ Joseph Fasel,^ John Middleditch,^ 
Mario Perez,^ Andrea Schmidt,^ John Singleton^ 

^Institute of Astronomy, University of Cambridge - UK 
^Clarendon Laboratory, Oxford University - UK 
^Los Alamos National Laboratory - USA 

Observational data imply the presence of superluminal electric currents in pulsar magnetospheres. 
Such sources are not inconsistent with special relativity; they have already been created in the 
laboratory. Here we describe the distinctive features of the radiation beam that is generated by a 
rotating superluminal source and show that 

(i) it consists of subbeams that are narrower the farther the observer is from the source: sub- 
beams whose intensities decay as 1/R instead of with distance (R), 

(ii) the fields of its subbeams are characterized by three concurrent polarization modes: two 
modes that are 'orthogonal' and a third mode whose position angle swings across the sub- 
beam bridging those of the other two, 

(iii) its overall beam consists of an incoherent superposition of such coherent subbeams and 
has an intensity profile that reflects the azimuthal distribution of the contributing part of 
the source (the part of the source that approaches the observer with the speed of light and 
zero acceleration), 

(iv) its spectrum (the superluminal counterpart of synchrotron spectrum) is broader than that of 
any other known emission and entails oscillations whose spacings and amplitudes respec- 
tively increase and decrease algebraically with increasing frequency, and 

(v) the degree of its mean polarization and the fraction of its linear polarization both increase 
with frequency beyond the frequency for which the observer falls within the Fresnel zone. 



We also compare these features with those of the radiation received from the Crab pulsar. 
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1. Introduction 

The rigid rotation of the overall distribution pattern of the pulsar emission reflects a radiation 
field E whose cylindrical components depend on the cylindrical coordinates (r, (p,z) and time t as 

Er^^^,{r,(p,z;t) = Er^^^,{r,(p - (Ot,z), (1.1) 

where co is the rotation frequency of the pulsar. Such a field can only arise from an electric current 
whose density j likewise depends on the azimuthal angle (p in the combination (p — cot only: 

jr,(p,z{f^(P,Z;t) = jr,(p,z{r,(p-0)t,z) (1.2) 

(see appendixes A and B of |[I]]). This property of the emitting current follows not only from 
the observational data, but also from the numerical models of the magnetospheric structure of an 
oblique rotator; it is found that any time-dependent structures in such models rapidly approach a 
steady state in the corotating frame [Q]. 

Unless there is no plasma outside the light cyUnder, therefore, the macroscopic distiribution of 
the emitting current in the magnetosphere of a pulsar should have a superluminally rotating pattern 
in r > c/ft) (where r is the radial distance from the axis of rotation and c is the speed of light 
in vacuo). Such a source is not inconsistent with special relativity. The superluminally moving 
pattern is created by the coordinated motion of aggregates of subluminally moving particles |^]. 
It has been experimentally verified, on the other hand, that such moving charged patterns act as 
sources of radiation in precisely the same way as any other moving sources of electromagnetic 
fields §110]. 

To illustrate the distinctive features of the emission from a superluminal source, we consider 
a polarization current whose distribution pattern rotates and oscillates at the same time: j = dV/dt 
for which 

Pr.(p^,{r,(p,z,t) = Sr.^^,{r,z)co^{m^))co^,{Q.t), -n<q><n, (1.3) 

and (p = (p — cot. Here, Pr,(p,z are the components of the polarization P in a cylindrical coordinate 
system based on the axis of rotation, s(r,z) is an arbitrary vector function with a finite support 
in r > c/ft), m is a positive integer, and Q. is an angular frequency whose value differs from an 
integral multiple of the rotation frequency co. This is a generic source: one can construct any 
distribution with a uniformly rotating pattern, P,-,(p,zif, ^,z), by the superposition over m of terms of 
the form Sr_(p_^{r,z,m) cos (m^). It also corresponds to laboratory-based sources that have been used 
in experimental demonstrations of some of the phenomena described below [Q, ^]. 

The results reported here are derived from the retarded solution of Maxwell's equations for the 
above current distribution (see [III ||, 0, 110, III], O]). 



2. The field generated by a constituent volume element of the source 

A superluminal source is necessarily volume-distributed [^]. However, its field can be built 
up from the superposition of the fields of its constituent volume elements which are point-like. 
Figure |l| shows that the waves generated by a constituent volume element of a rotating superluminal 
source possess a cusped envelope and that, inside the envelope, three wave fronts pass through any 
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Figure 1: Left: envelope of the spherical wave fronts emanating from a superluminally moving source 
element (S) in circular motion. The heavy curves show the cross section of the envelope with the plane of 
the orbit of the source. The larger of the two dotted circles designates the orbit and the smaller the light 
cylinder Middle and right: the spacetime diagrams showing the intersection of the trajectory of the source 
point S with the past light cone of the observation point P when P lies inside, and on the cusp of, the envelope 
of wave fronts, respectively. 





(a) (b) 

Figure 2: Three dimensional views of the envelope (a) and its cusp (b). 



given observation point simultaneously. This reflects the fact that the field inside the envelope 
receives simultaneous contributions from three distinct values of the retarded time (see Fig. [l]). On 
the cusp of the envelope, where the spacetime trajectory of the source is tangent to the past light 
cone of the observer (Fig. all three contributions toward the value of the field coalesce [||, ^, p^ ]. 

On this cusp (caustic), the source approaches the observer with the speed of light and zero 
acceleration at the retarded time, i.e. dR{t)/dt = —c and d^R{t)/dt^ = 0, where R{t) = \x{t) — xp\ 
is the distance between the source point x{t) and the observation point xp. As a result, the interval 
of emission time for the signal carried by the cusp is much longer than the interval of its reception 
time [|]. 

A three-dimensional view of the envelope of wave fronts and its cusp is shown in Fig. ^. The 
two sheets of the envelope, and the cusp along which these two sheets meet tangentially, spiral out- 
ward into the far zone. In the far zone, the cusp lies on the double cone dp = arcsin[c/(rG))], dp = 
71 — arcsin[c/(rft))], where {Rp,dp,(pp) denote the spherical polar coordinates of the observation 
point P. Thus, a stationary observer in the polar interval arcsin[c/(rft))] < dp < n — arcsin[c/(rft))] 



receives recurring pulses as the envelope rotates past him/her [ 1 1 1. 

Figure ^ shows the radiation field generated by the rotating source element 5 on a cone close 
to the cusp, just outside the envelope. Not only does the spiralUng cusp embody a recurring pulse, 
but the plane of polarization of the radiation swings across the pulse [jl^]. 

As a cosequence of the multivaluedness of the retarded time, three images of the source are 
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Figure 3: Polarization position angles and field strengths on the cone Op = n/12 outside the envelope for a 
source with rco = 2. 




Figure 4: An observer P inside the envelope detects three images (/i, h, I3) of the source S simultaneously. 



observeable inside the envelope at any given observation time (Fig. ^). The waves that were emitted 
when the source was at the retarded positions /i , I2 and /a in Fig. ^ are all received simultaneously 



at the observation point P. These images are detected as distinct components of the radiation [13]. 



Field strengths and polarization position angles of the three images (radiation modes) close to 
the cusp are shown in Fig. Two modes dominate everywhere except in the middle of the pulse. 
Moreover, the position angles of two of the modes are 'orthogonal' and that of the third swings 
across the pulse bridging the other two. The constructive interference of the emitted waves on the 
envelope (where two of the contributing retarded times coalesce) and on its cusp (where all three 
of the contributing retarded times coalesce) gives rise to the divergence of the field of a point-like 
source on these loci M\. Here we plot the spatial distribution of the field excising the narrow regions 



in which the magnitude of the field exceeds a certain threshold [12] 
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Figure 5: Left; the relative strengths of the three radiation modes as observed near the cusp on a sphere 
of large radius. The total field strength (black) and strengths of the underlying contributions from the three 
images of the source (green, red, blue) are shown for a source with rco = 1.1 and an observation point that 
sweeps a small arc of the circle RpCo/c = lO'", 9p ~ Tt/l.l, crossing the envelope near the cusp. Right: the 
corresponding position angles of the contributions from the three retarded times (green, red, blue) are shown 
relative to one another and to that of the total field (black); the position angles of the dominant contributions 
are shown with open triangles, and those of the weakest contributions with filled triangles. 



3. The nondiffracting subbeams comprising the overall beam 

The dominant contribution towards the field of an extended source comes from a thin filamen- 
tary part of the source that approaches the observer, along the radiation direction, with the speed 



of light and zero acceleration at the retarded time []11|]. For an observation point P in the far zone 
with the coordinates {Rp, dp, (pp), this filament is located at r = (c/ft)) esc dp, cp = (pp + 37r/2 and 
is essentially parallel to the rotation axis (Fig. ^). The collection of cusps of the envelopes of wave 
fronts that emanate from various volume elements of the contributing filament form a subbeam 
whose polar width is nondiffracting: the linear dimension of this bundle of cusps in the direction 
parallel to the rotation axis remains the same at all distances from the source, so that the polar angle 
5 dp subtended by the subbeam decreases as Rp^^ with increasing Rp (see Fig. ^. 

In that it consists of caustics and so is constantly dispersed and reconstructed out of other 
waves, the subbeam in question radically differs from a conventional radiation beam [j8|]. The 
narrowing of its polar width (as Rp^^) is accompanied by a more slowly diminishing intensity (an 
intensity that diminishes as Rp^ instead of Rp^ with distance), so that the flux of energy across 



its cross sectional area remains the same for all Rp [|1 1|]. This slower rate of decay of the emission 
from a superluminally rotating source has been confirmed experimentally [§, ^. 

The contributing part of an extended source (the filament that approaches the observation point 
with the speed of light and zero acceleration) changes as the source rotates (see Fig. |6|). In the case 
of a turbulent plasma with a superluminally rotating macroscopic distribution, therefore, the overall 
beam within which the narrow, nonspherically decaying radiation is detectable would consist of an 
incoherent superposition of coherent, nondiffracting subbeams with widely differing amplitudes 
and phases (similar to the train of giant pulses received from the Crab pulsar [|14|]). 

The overall beam occupies a solid angle whose polar and azimuthal extents are independent 
of the distance Rp. It is detectable within the polar interval arccos[(r<G)/c)^'] < \dp — n/2\ < 
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Figure 6: Schematic illustration of the light cylinder, the filamentary part of the source that approaches the 
observation point with the speed of light and zero acceleration at the retarded time, the orbit r = c/{(Osin dp) 
of this filamentary source, and the subbeam formed by the bundle of cusps that emanate from the constituent 
volume elements of this filament. 



arccos[(r>G)/c)^'], where [r<,r>] denotes the radial extent of the superluminal part of the source. 
The azimuthal profile of this overall beam reflects the distribution of the source density around the 
cylinder r = c/(G)sin dp), from which the dominant contribution to the radiation arises 1 11]. 



4. Superluminal counterpart of synchrotron spectrum 

The spectrum of the radiation emitted by a rotating superluminal source is oscillatory with 
oscillations whose spacing increase with frequency ^. While the Bessel function describing 
synchrotron radiation has an argument smaller than its order and so decays exponentially with 
increasing frequency, the Bessel function encountered in [^], whose argument exceeds its order, 
is an oscillatory function of frequency with an amplitude that decays only algebraically. Figure ^ 
shows that the spacing of the emission bands in the spectrum of the Crab pulsar fit the predicted 
oscillations for an appropriate choice of the single parameter Q./(0. The value of this parameter, 
thus implied by the data of [|^], places the last peak of the oscillating spectrum at a frequency 
(~ Q? /ojP') that agrees with the position of the ultraviolet peak in the spectrum of the Crab pulsar. 



By inferring the remaining adjustable parameter m in Eq. ( |1.3| ) from the observational data and by 
mildly restricting certain local properties of the source density s, we are thus able to account for 
the continuum spectrum of the Crab pulsar over 16 orders of magnitude of frequency (Fig. |8|). 

At radiation frequencies higher than ~ 10^^ Hz, the Earth falls within the Fresnel zone of the 
Crab pulsar. For an observer in this zone, the emission arises from a narrower radial extent of the 
source, and so has a higher degree of mean polarization, the higher the frequency. The degree of 
circular polarization of such a high-frequency emission decreases with increasing frequency v as 
v^^/^ , so that this emission is essentially 100 per cent linearly polarized at all phases, including 
off-pulse phases (see Eq. (15) of [|I|]). Note that, in this model, the intensity and polarization of the 
off-pulse emission, too, reflect the distributions of density and orientation of the emitting current 



around the cylinder r = c/(g) sin Qp) within the pulsar magnetosphere [[111]. 

Spectral differences of the different parts of the source that contribute to the emission detected 
at different rotation phases can result in the frequency dependence of the position angle of the 
high-frequency radiation. Thus, one can have a linearly polarized high-frequency emission whose 
position angle is independent of phase but dependent on frequency (as observed in the Crab pul- 
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Figure 7: The predicted oscillations of the spectrum of the emission for co/(2n) = 30 Hz and Q./co = 
1.9 X 10"* (left) have the same spacing as those of the emission bands in the observed spectrum of the Crab 
pulsai- (right). 




logio 1^ (Hz) 



Figure 8: The data points show the continuum spectrum of the Crab pulsar [|T6|. In the model, the recovery 
of intensity at the ultraviolet peak is caused by resonance with mco/ {2n) ~ 3 THz. The spectral break at 
~ 10^^ Hz reflects a transition across the boundary of the Fresnel zone (see [|lj). 



sar [17]): independent of phase because the emitting current density can have a single non-zero 
component (e.g. the component parallel to the roation axis) at all points around the cylinder 
r = c/{cosmdp), but frequency dependent because the distribution of the source as a function of (p 
can depend on frequency 



5. Comparison with observations 

Characteristic features of the radiation generated by the superluminal emission mechanism 
described above suggest 

(i) that the extreme values of the brightness temperature (~ 10^' K), temporal width (~ 1 ns), 
and source dimension (~ 1 m) of the giant pulses received from the Crab pulsar arise from 
the nonspherical decay of the intensity of this radiation with distancce, 
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(ii) that the observed micro- and nanostructures of the radiation beams received from the Crab 
pulsar arise from an incoherent superposition of sets of coherent, nondiffracting subbeams, 

(iii) that the breadth of the spectrum of the radiation received from the Crab pulsar, and its propor- 
tionately spaced emission bands, reflect the algebraic (instead of exponential) rate of decay 
and oscillations of the superluminal counterpart of synchrotron spectrum, 

(iv) that concurrent 'orthogonal' polarization modes with swinging position angles arise from 
the multiple images of a compact source and so should be associated with individual giant 
pulses, 

(v) that the unpulsed part of the emission from the Crab is also generated by the pulsar itself, 
and 

(vi) that the high degree of polarization of the high frequency radiation from the Crab pulsar 
reflects the decrease in size of the contributing part of the source with increasing frequency 
and a corresponding increase in degrees both of mean polarization and of Unear polarization 
with increasing frequency. 
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